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Persistence of a Non-Equilibrium State: Observation of a Boltzmannian Special Case 
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Well before the atomistic nature of matter 
was experimentally established, Ludwig Boltz- 
mann’s audacious effort to explain the macro- 
scopic World of human experience in terms of 
the workings of an unseen microscopic world met 
with vigorous oppositiorP^. A contentious point 
was the problem of irreversibility: the micro¬ 
scopic equations of motion are reversible, yet fric- 
tion and viscosity cause things always to slow 
down and warm up, never to speed up and cool 
dowrP. What was worse, Boltzmann himself dis- 
covered that his transport equation predicts spe¬ 
cial cases in which gases never come to thermal 
equilibrium, a particular example being that the 
monopole “breathe” mode of gas will never damp 
if it is confined in 3D to a perfectly isotropic 
harmonic potentiaP. Such absences of damp- 
ing were not observed in nature. Nondamping 
of a monopole mode in lower dimensional Sys¬ 
tems has only very recently been observed, us- 
ing cold atoms. Kinoshita et ah® and Chevy 
et ah® have experimentally observed suppressed 
relaxation in highly elongated geometries. The 
difflculty in generating sufficiently spherical har¬ 
monic confinement for ultracold atoms, however, 
has meant that Boltzmann’s fully 3D, isotropic 
case has never been observed. With the develop- 
ment of a new magnetic trapP capable of produc- 
ing near-spherical harmonic confinement for ul¬ 
tracold atoms, we have been able to make the first 
observation of this historically significant oddity. 
We observe a monopole mode for which the col- 
lisional contribution to damping vanishes, a long- 
delayed vindication for Boltzmann’s microscopic 
theory. 

The Boltzmann equation determines how the phase- 
space distribution of a gas, / (r, v,t), evolves as a func- 
tion of binary collisions between particles with mass, m, 
in the presence of an external force, F 

f + ^ • + - • Vv/ = IcoU [/] . (1) 

at at m 

The collision integral, Icoiu describes how populations at 
the same location, r, but differing velocities, v, redis- 
tribute to two new velocities, v'. These local, pairwise 
collisions conserve momentu m a nd energy. For the ex- 
plicit form of the integral, se^^. 

The collision integral vanishes whenever the product 
of two single-particle distributions is identical directly 
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before and after a collision, in other words, when 

/(Vl)/(V2) = /(Vi')/(V2')- (2) 

This equality typically implies that the phase-space dis¬ 
tribution is time-invariant and the gas has reached equi¬ 
librium. However, when energy, momentum, and total 
atom number are conserved, equation is generically 
satisfied by distributions of the form 

/ (r, V, t) cx C/(r)e ^ ('3) 

where U (r) contains information about the external con- 
fining potential, kg is the Boltzmann constant, the tem¬ 
perature T(t) is time-dependent and r](r,t) is an arbi- 
trary function of space and time. While these distri¬ 
butions, known as “local equilibrium distributions”, al¬ 
ways cause the collision integral to vanish, they in gen- 
eral do not satisfy equation By constraining the 

local equilibrium distribution so that ? 7 (r,t) = 0 and 
dT/dt = 0, the distribution becomes a valid solution of 
equation 0 known as the Maxwell-Boltzmann distribu¬ 
tion: 

/ (r, v) (x U(r)e~ . (4) 

But certain potentials exist where equation Q is satisfied 
by non-equilibrium distributions, in which case the time- 
dependence in equation 0 remain^Sliol Qj^g gf t^gge 
cases is the 3D isotropic harmonic potential with a solu¬ 
tion corresponding to a spherically symmetric “monopole 
mode”, where temperature and cloud size oscillate, with 
opposite phase, in tim^i^HIIl Because the temperature is 
oscillating in the absence of heat conduction, it is conve¬ 
nient to call the temperature in equation ([^ a “kinetic 
temperature”, Tk (t) , and deline a “spatial temperature”, 
Ts(t), which determines the variation in cloud size. The 
average temperature (Ifc(t) -|- Ts{t))/2 is constant; the 
breathing dynamics are analogous to the oscillatory ex- 
change between kinetic and potential energy that occurs 
in simple harmonic motion. While the structure of equa¬ 
tion Q implies a quasicontinuous distribution / and thus 
very large atom number N, we show in the methods sec- 
tion that the resuit of vanishing damping is preserved as 
N increases from 1 to 2, and onto an arbitrary meso- or 
macroscopic number. 

This strange absence of damping holds for the 
monopole mode, but not necessarily for other collective 
modes. For the quadrupole mode, a mode in which the 
radial and axial widths oscillate 180° out of phase, cross- 
dimensional coupling from collisions causes damping. In 
the limit of an interatomic collision rate, ^coU^ that is 
much smaller than the trap frequency, the quadrupole 
mode in an isotropic harmonic potential is predicted to 
damp at a rat^^^^ 
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We will measure quadrupole damping rates as a baseline 
for comparison with measured monopole damping rates. 

The experiment is per formed with ®^Rb atoms evap- 
oratively cooled in a magnetic trap with har- 

monic confinement at frequency w = 27r(9.03(2) Hz) and 
equipped with additional magnetic coils that permit the 
six distinet parameters of a 3D quadratic potential to 
be adjusted independentljl^. We measure dipole slosh- 
ing motion of atoms in the trap to determine i, j, k, the 
principle axes of the confining potential, and their asso- 
ciated trapping frequencies, tOi, ujj, ujk- We characterize 
the residual asphericity, {ujmax - where 

Wmm and <D are respectively the maximum, minimum and 
mean of Wi, Wj, The residual asphericity drifts with 
time so we periodically retune and recharacterize the trap 
to keep asphericity small, typically less than 0.002. To 
minimize the undesireable mean-held potential, we work 
at temperature T well above the Bose-Einstein conden¬ 
sate transition temperature, Tc, between 2Tc <T < 3Tc. 

We selectively drive monopole (quadrupole) motion by 
symmetrically (asymmetrically) modulating the strength 
of the confinement about its mean value. The cloud is 
then allowed to evolve freely in the spherical trap be- 
fore it is n on-de structively imaged using phase-contrast 
microscopjff^E^. For each cycle of the experiment, six 
images are taken of the cloud along two orthogonal axes 
with an interval of 17 ms in order to sample roughly 1.5 
oscillation periods. Cloud widths along each dimension, 
o'i,j,k, are determined using Gaussian surface fits of in- 
dividual images in order to determine the amplitude of 
the monopole and quadrupole modes. Amplitudes of the 
monopole and quadrupole distortion are scaled by the 
average width of the cloud during one cycle, given by the 
following relations 


af + a'j + al 

A _ ^ J K 

■^M / 2 I 2 ’i 2\ ^ 
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Aq = 


2al -af-a] 

+^r+4y 


( 7 ) 


Oscillation amplitudes are determined by fitting a cycle 
of oscillation in cloud width, obtained from each experi- 
mental run, with a fixed frequency sine wave as indicated 
by the solid lines in figure 

Suppressed damping of the monopole mode relative to 
the quadrupole mode can be seen in the sample data 
in figure and although the damping rate is small, it 
is nonzero. We characterize the monopole damping by 
comparing with quadrupole mode damping rates, which 
are expected to vary linearly with collision rate. By ad- 
justing the evaporation parameters in our experiment, 
we can tune N, T, and the collision rate of the sam¬ 
ple, and then alternately drive quadrupole or monopole 
modes. A direct comparison of quadrupole and monopole 
damping rates in a near-spherical trap is shown in hgure 
The dependence of quadrupole damping on collision 



0 200 400 600 800 1000 1200 1400 1600 

Time (ms) 


FIG. 1. Sample data for a driven quadrupole mode and 
monopole mode in a spherical trap with residual asphericity 
less than 0.002 and a collision rate of 7.4(3) s“^. Solid lines 
on the quadrupole data indicate a typical fitting procedure 
where individual periods taken in a single run are fit with an 
undamped sine wave to extract an instantaneous amplitude. 
The oscillation amplitudes at various cycles are in turn fit to 
an exponential decay to extract the damping rates shown in 
figure Only a small subset of the monopole mode data is 
shown and the full set spans 30 seconds. Random observ- 
able scatter in these points is predominantly due to small, 
irreproducible fluctuations in initial conditions. 

rate is Tq = (4.9(1))“^7coH, which is in good agreement 
with equation ([^. The small amount of residual damp¬ 
ing in the monopole mode is independent of collision rate 
and much smaller than the damping in the quadrupole 
mode, as expected. This, then, is the special-case ex- 
ception that proves the general rule of damping in the 
Boltzmann equation. 


g 

'a 


Monopole ■ 
Quadrupole • 


y 


.4' 

■ 

yf 


- 



- 

■ ■■ 

■ 

■ 


0 2 4 6 8 10 

Collision rate (s“^) 


FIG. 2. Monopole and quadrupole damping rates as a 
function of interatomic collision rate in a near-spherical trap. 
The dotted line indicates the predicted quadrupole damping 
rates given in equation (§. Residual asphericity was typically 
around 0.001, but always less than 0.0026. Cloud full width 
at half maximum was approximately 115 fim. Error bars in 
the monopole data are smaller than the data points. There is 
a multiplicative uncertainty in the collision rate of 10%. 
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In order to understand the source of residual monopole 
damping, we note that Boltzniann’s resuit hinges on the 
assumption that the potential is both isotropic and har- 
monic. An actual physical System can never satisfy both 
of these conditions perfectly, and the remainder of this 
letter is devoted to a discussion of the effects that small 
anisotropies and anharmonicities have on the monopole 
damping rate. 

Certain subtleties arise for gases in the collisionless 
limit when anisotropies in the potential are small enough 
that trap frequencies differ by less than a few percent. In 
a totally collisionless System, oscillations along the Prin¬ 
cipal axes of the trap are fully decoupled and monopole- 
or quadrupole-like oscillations are undamped. If the Prin¬ 
cipal trap frequencies differ such that = uij ^ de- 
phasing occurs between oscillations along different prin- 
cipal axes and energy exchange between pure monopole 
motion and pure quadrupole motion occurs with a period 
given by 


Tmq = 


TT 

|wfc - Wil ■ 


( 8 ) 


When collisions are included, the two modes become cou- 
pled and, as the population in the quadrupole mode in- 
creases, so does the dampinjl^. This effect can be seen 
when Tmq < l/bg, where multiple oscillations between 
monopole and quadrupole modes occur. Data in figure 
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FIG. 3. Monopole and quadrupole oscillation amplitudes 
in an anisotropic trap where the residual asphericity is ap- 
proximately 0.02. The data points show the amplitudes of 
individual oscillations of the instantaneous distortion, such as 
those indicated by the black lines in Fig.j^ A monopole mode 
is initially driven and oscillations between the monopole and 
quadrupole modes at a frequency Au can be seen. The dotted 
line fitting functions are (cos^(Aut-|-<5) + C). Ampli¬ 
tudes are determined from sine-wave fits to individual oscilla¬ 
tion measurements obtained from a single experimental cycle. 
The error bars are calculated from scatter in the sine-wave fits 
of data taken at repeated evolution times. 


show oscillations between monopole and quadrupole 
modes resulting from an initial monopole drive in a trap 
with residual asphericity of approximately 0.02. En¬ 
ergy transfers back and forth between the individual 
modes and the damping rates for both modes are nearly 
equal, with a mean value of T = 0.36(4) s“^. The 
collision rate is roughly 3.7 s“^ leading to an expected 
quadrupole damping rate of 0.74 s“^ in a spherical trap, 
which is twice the value of the measured damping rate 
in the anisotropic trap. This is no surprise because the 
quadrupole mode is effectively populated only half of the 
time, leading to the factor of 2 decrease in the damp¬ 
ing rate. If we decrease the amount of anisotropy such 
that Tmq >> l/bg, the energy in the quadrupole mode 
damps before it can fully couple back into the monopole 
mode. One can see this effect in the data for the very 
spherical case shown in figure where the quadrupole 
mode damps before it can exchange with the monopole 
mode. The data in figure were taken in traps with 
residual asphericities ranging from 0.0005-0.0026, corre- 
sponding to 10.7 s < Tmq/‘2. < 55.6 s. The typical re- 
laxation time for the monopole mode in these traps is 
tm = 7(2) s with no systematic dependence on the resid¬ 
ual anisotropy observed. Thus, some other physical effect 
must provide the dominant source of residual monopole 
damping. 

We now come to the second condition of Boltzmann’s 
resuit, which requires that the potential be harmonic, 
and discuss the effect of anharmonic perturbatione to our 
trapping potential as a source of damping. Amplitude- 
dependent frequency shifts caused by anharmonic pertur- 
bations lead to dephasing of individual particle trajec- 
tories, effectively damping the collective monopole am- 
plitude. Moreover, the anharmonic correctione to our 



FIG. 4. Damping of the monopole mode in a near-spherical 
trap as a function of the spatial extent of the atom cloud. 
The fractional amplitude of the excitation is the same for all 
points. Consistent with anharmonicity-induced damping, we 
see the observed damping rate decreases rapidly for smaller 
clouds, but the trend of the data suggest that even for the 
smallest clouds (such as those used in figure the residual 
damping observed may be due to the onset of anharmonic ef¬ 
fects. The error bars are calculated using the same method 
described in Fig.[^and the error in the cloud width measure¬ 
ments were typically 2%. 
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potential are asymmetric, giving rise to an amplitude- 
dependent anisotropy. A calculation of the expected 
damping rate that takes into account all of the relevant 
anharmonic corrections is difficult. But the effect can 
be explored experimentally by measuring the monopole 
damping as a function of cloud size as shown in figure 
The first point in figure represents the average of the 
cloud size and monopole damping data in figure The 
trend of the data in hgure suggests that the residual 
damping seen even at the smallest cloud size of 115 ^m 
is already due to the onset of anharmonic effects. Unfor- 
tunately, we are unable to work with smaller, and thus 
colder, clouds due to limitations in the signal-to-noise ra¬ 
tio of our imaging system, and the need to keep T > 2Tc. 

In this paper, we present an experimental verfication 
of the absence of damping for the monopole mode of a 
thermal gas in an isotropic harmonic potential. While 
the damping is highly suppressed, the small but finite re- 
laxation of the monopole mode is an artifact of small an¬ 
harmonic perturbations to our trap, which decrease with 
cloud size. We find that, in the limit of zero anharmonic 
shifts, the damping of the monopole mode vanishes, as 
predicted by Boltzmann in 187fP. 


METHODS 


The undamped nature of the monopole mode is found 
by calculating the evolution of the square radius of the 
cloud and can be derived in various way^^. It is instruc- 
tive to see how the monopole nondamping resuit can be 
built-up starting from = 1, 2, etc... In spherical 
symmetry, the radial motion of a single particle of mass 
TO, energy E, and angular momentum L is governed by 
the effective potentiaP^ 


TA ^ 1 2 2 

so that the radial force is 


( 9 ) 




( 10 ) 


and the kinetic energy is 




( 11 ) 


We note that (Pr‘^/dt‘^ = 2{dr/dt)"^ + 2r(Pr/dt^, and sub- 
stituting ([To| and pll yields the differential equation for 


Ar» =-nqr-- rS , (12) 

where fl = 2 uj and Tq = Ej{muj'^). So the square radius 
undergoes sinusoidal oscillations, or “monopole breathe”, 


around its mean value Cg at a frequency of 2a;. If there 
are two particles, 1 and 2, each with individual values of 
E, L, and r^, each particle will oscillate at 2uj. Taking the 
sum of their respective differential equations (121 yields 


dt'^ 


{r! - rl) 


(13) 


where their combined square radius, -|- r|, oscil- 

lates around its mean value, Cg^ = {Ei + E 2 )/. The 
magnitude of the collective breathe motion depends on 
the magnitude and relative phase of the individual parti¬ 
cle trajectories. These individual quantities will abruptly 
change in the event of a collision. Assuming the collisions 
are local, ri, r 2 , and thus will not change from the in¬ 
stant before to the instant after the collision. Similarly, 
momentum and energy conservation imply that and 
Tgj are unchanged by the collision. 

These three continuities imply that the parameters and 
boundary conditions of equation (131 are matched di- 
rectly before and after a collision. This ensures that 
neither the magnitude nor phase of the oscillation will 
change as the resuit of a pairwise collision. If we instead 
consider N atoms where 


N 



i=l 


(14a) 
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(14b) 


One can see that the monopole mode is left 
unperturbed—and in particular undamped —by local, 
pairwise, momentum-, energy-, and number-conserving 
collisions. This argument is robust to quantum 
statistics—Bose or Fermi—and, interestingly, in Ref. HHl 
it is shown that a l/r^ term in the potential also pre- 
serves monopole motion. 

The above approach is consistent with a system in 
the hydrodynamic limit, with total number of atoms so 
large that the function /(r,v,t) is essentially continu- 
ous. However, in the hydrodynamic limit, mean-field 
effects can come into play in which case the monopole 
frequency is shiftecP^. In our experiment, the total num¬ 
ber of atoms, N, is only a few hundred thousand, and 
the mean-free path is large compared to the spatial ex- 
tent of the sample—we are not in the hydrodynamic limit 
in any sense of the word—and mean-field effects can be 
neglected. 
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